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We investigate the properties of mesons with the exotic = 1“'*’ quantum numbers. Starting 
out from the light-quark domain, where the tti states are used as references, we predict the masses 
of analogous quarkonia for cc and hb configurations. We employ a covariant Dyson-Schwinger-Bethe- 
Salpeter-equation approach with a rainbow-ladder truncated model of quantum chromodynamics. 
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I. INTRODUCTION 

After its origin several decades ago, hadron spec¬ 
troscopy has seen increased interest in recent years. In 
particular, states not predicted or readily explained by the 
comprehensive quark-model calculations of the past, e.g., 
[T] have sparked both theoretical investigations as well 
as experimental searches at modern and future hadron- 
physics facilities; foremost, in the present manuscript we 
focus on mesons with exotic quantum numbers, which 
have been reviewed in great detail, e. g., very recently [5]. 

The approach we chose is the coupled Dyson-Schwinger- 
Bethe-Salpeter-equation (DSBSE) framework for the in¬ 
vestigation of quantum chromodynamics (QCD), a well- 
established covariant variant [m from a set of continuum 
quantum-field-theoretical methods [MD] complementary 
to the well-known lattice-regularized QCD [IIHS]. 

Concretely, our present work is built on top of recent 
studies of heavy quarkonia mui], to which we add an 
immediate focus on the exotic meson quantum numbers 
= 1 thus approaching the quarkonium spectrum 
from a slightly different perspective. In particular, the 
present study reaches out to the light-quark domain for 
anchoring in addition to the conventional quarkonium 
spectra in order to make use of experimental evidence 
presently at hand for the tti states. Our setup employs 
the basic but symmetry-preserving rainbow-ladder (RL) 
truncation to study mesons by solving the quark Dyson- 
Schwinger equation (DSE) coupled to the meson qq Bethe- 
Salpeter equation (BSE). 


II. SETUP AND MODEL 

While we perform an RL-truncated study of Landau- 
gauge QCD in Euclidean space using the homogeneous 
BSE, we refer the reader to alternate routes of studies such 
as ours via a brief list of exemplary references to Coulomb- 
gauge QCD [IMT^ . the BSE in Minkowski space [201122] . 
an analogous but more general inhomogeneous vertex 
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BSE and studies beyond RL truncation 

from which the relevant literature can be grasped. 

The strength of the setup chosen herein is its com¬ 
bination of computational and conceptional feasibility 
I3I1ISH with the important features inherent to a system¬ 
atic setup of a nonperturbative quantum-field-theoretical 
method. In particular, in QCD the implementation of 
chiral symmetry and its dynamical breaking as well as a 
realization of confinement and meson properties expected 
in the heavy-quark limit are paramount. Our study meets 
these requirements via properties demonstrated of RL 
truncated DSBSE studies with dressed quarks before, 
namely satisfaction of the axial-vector Ward-Takahashi 
identity [3M3^ , quark-confinement |36l [37] , and adequate 
behavior towards the heavy-quark limit [551140] . 

We use the homogeneous qq BSE in RL truncation 
which reads 

r(p; P) = g{{p-qf) Dl,,{p-q) 7 ^ x(g; P) 7 ^ 
xiq;P) = S{q+)Tiq-,P)S{q.), (1) 

where q and P are the quark-antiquark relative and total 
momenta, respectively, and the (anti)quark momenta are 
chosen as q± = q± P/2. The renormalized dressed quark 
propagator S{p) is obtained from the corresponding DSE 

= (^2*7 ■ P + ZilTlqip.)) + S(p) , 

S(p) = dU(p-q) In S{q) lu- (2) 

E is the quark self-energy, rriq is the current-quark mass 
chosen at a renormalization scale /i, represents the 
free gluon propagator and Xu is the bare quark-gluon 
vertex’s Dirac structure. Dirac and flavor indices are 
omitted for brevity, = f d'^q/ifli^Y denotes a trans- 
lationally invariant regularization of the integral, with the 
regularization scale A sa¬ 
in a given truncation the effective interaction Q needs 
to be specified according to the aims of the study. We use 
the well-investigated and phenomenologically successful 
parameterization of Ref. [55] , which reads 

g(s) ^ 47r^-P -sjuP , 47 r T{s) 

s « ^l/ 21 n[r+(l+s/A^cD)^]’ ^ ^ 
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The parameter oj corresponds to an effective inverse 
range of the interaction, while D acts like an overall 
strength of the first term; they determine the intermediate- 
momentum part of the interaction, while the second term 
is relevant for large momenta and produces the correct 
one-loop perturbative QCD limit. We note that J^{s) = 
[1 — exp(— s/[4toj])]/s where rrit = 0.5 GeV, r = e^ — 1, 
N{ = 4, = 0.234 GeV, and 7 ^ = 12/(33 - 2Vf), 

which is left unchanged from Ref. [42]. 

The parameters considered independent and free in 
our study are ui and D (in addition to the current-quark 
masses), which we already used previously to achieve a 
surprisingly good phenomenological description of meson 
spectra mills] in the heavy-quark domain. In the present 
study we refit the pair oj and D using the same quark 
masses as in mi EH in order to specifically target the 
= 1 channel. Further details are specified below. 


III. EXOTIC QUANTUM NUMBERS 

Exotic quantum numbers have been addressed within 
the DSBSE approach several times before mi 13H1 ElHin] , 
but have been in the actual focus of the work only in 

mm?]. 

In a short critique of Ref. m we would like to mention 
here in particular the claim that uj and D are not indepen¬ 
dent so that “one can expect computed observables to be 
practically insensitive to uj on the domain uj G [0.4, 0.6] 
GeV.” In our opinion, this is an unnecessary and over- 
restrictive assumption which actually prevents part of the 
possibilities a complete set of results may have. We note 
that there is a difference between having a one-parameter 
model at hand like in older studies to investigate effects 
the long-range behavior of the strong interaction might 
have on states of various kinds [Ml IMII53] on one hand, 
or attempting a complete study of the possibilities within 
the parameter space of a given model interaction on the 
other hand. While in m the former is present and, in 
fact, a different value for the product of ujD is used for 
excited and exotic states to attempt a description along 
the lines of beyond-RL premises, the latter is not, and so 
any conclusions drawn there are necessarily incomplete. 
For example, general statements about excited and exotic 
states are made in m based on light-quark results only, 
while heavy mesons are not considered, nor are those with 
axial-vector and tensor quantum numbers. In our present 
work, on the contrary, we show that axial-vector states 
may be the key to a proper understanding of, e. g., the 
exotic vector states in the light-meson spectrum. 

In another short critique, namely of Ref. [45] . we remark 
that, after an in-depth discussion of the DSBSE setup and 
it’s capabilities to describe states with exotic quantum 
numbers, the authors employ a strategy similar to ours. 
In particular they focus on states with J = 1, together 
with the obligatory pseudoscalar ground state, and obtain 
reasonable agreement with experimental data for vector, 
axial vector, and exotic vector ground states. Since the 


choice of a separable kernel in the BSE limits the resulting 
equations and bound-state amplitudes, this set of results 
cannot easily be generalized or extended. For example, 
only two solutions are supported in the 1 channel 
which also appear rather close together in mass, and an 
immediate comparison to the RL setup employed herein 
is difficult. At the time, possible interpretations as to 
why the results in |45j compared more favorably to the 
experimental situation than RL studies were mainly that 
there was lesser, and thus more appropriate, strength 
in the separable as compared to the RL kernel for light 
quarks. 

With these two DSBSE studies of states with exotic 
quantum numbers in mind, we now very briefly sketch a 
few relevant issues before outlining our strategy towards 
meaningful statements in our approach. 

In the DSBSE approach, states with exotic quantum 
numbers appear naturally in the setup of the quark- 
antiquark BSE. Owing to the ubiquity of the scalar prod¬ 
uct of the quark-antiquark relative and total momenta, 
q • P, and its negative C-parity, the appearance of such 
states is a direct result of the dynamics and interaction 
considered. While in lattice QCD or constituent-quark 
models some explicit inclusion of gluonic degrees of free¬ 
dom leads to the exotic value of G, a similar mechanism 
includes gluonic dynamics here implicitly] for more de¬ 
tails and a thorough discussion, see [45] . For the present 
purpose and possible interpretations of our results, it is 
mainly important to remark here that our setup does not 
actually imply a non-hybrid interpretation. 

Regarding the experimental situation of states seen in 
the 1 channel, we are aware of an ongoing discussion 
about the nature and properties of the tti states as they 
are listed in |54j . However, improvements are expected 
due to the wide efforts undertaken both experimentally 
and theoretically to make the overall quality and under¬ 
standing of the hadron spectrum even better [^. For the 
purposes of our argument, we assume enough credibility 
of experimental information on the 7ri(1400) and 7ri(1600) 
states to warrant our hypothesis, which is certainly the 
case. The 7ri(2015) is neither needed for our fit setup, nor 
matched by any of our resulting states (our results are too 
high to match it) and so is omitted from our discussion. 


IV. STRATEGY 

Naturally, a study using a truncated set of equations 
like ours has to be taken with a grain of salt. As a result, 
we employ the following strategy in pursuit of meaningful 
predictions for 1 ^ meson masses in the quarkonium 

spectrum: While an overall successful description of states 
with J < 2 in RL truncation has been proven feasible in 
the heavy-quark domain there is no reason a priori 
to assume similar success in the light-quark sector. On the 
contrary, a comprehensively successful phenomenological 
description of the same quality, even in the isovector case, 
would be a surprise. Nonetheless, subsets of observables 
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FIG. 1. Light isovector quarkonium spectra of mesons with J < 2 including exotic quantum numbers, fitted to exotic-vector- 
related splittings. Calculated values (blue circles) are compared to experimental data (red boxes) from |54| . 


can be targeted by informed and suitably eclectic studies, 
which is our tactic herein as well. 

In particular, as already indicated above, we make use of 
the experimental information available on the two lowest- 
lying exotic TTi states, the 7ri(1400) and 7ri(1600), to test 
and constrain those meson mass splittings that may be ap¬ 
propriate, i. e., well-correlated to identify parameter sets 
potentially successful in the description of the 1 channel 
also for the other quark masses. To find these, we use com¬ 
parisons of calculated and experimental mass splittings 
and of sets of these on a grid oj S {0.4,0.5,0.6, 0.7,0.8} 
GeV X D G 10.9,1.1,1.2,1.3,1.4,1.5,1.7} GeV^ in the 
isovector case. The pion mass is always fixed to its ex¬ 
perimental value, which sets the current-quark mass in 
Eq. to rriu/d = 0.003 GeV (given at a renormalization 
poinf/j, = 19 GeV). We calculated meson masses for all 
quantum number sets for J < 2 on the grid and 
investigated the corresponding for a number of com¬ 
parisons. A more detailed account of this strategy as well 
as an illustration is presented in |43j . 

To confirm that our region of interest is indeed well- 
captured by the grid defined above, we investigated the 
TTi states and found that a combination of the radial 
splitting between the ground and first excited 1“+ states, 
together with the two splittings of each 1 ^ state to the 
ground-state pseudoscalar meson (pion) indeed provides 
a region of low in the interior of our grid, which is the 
necessary prerequisite. We note that the best fit according 
to this configuration is given by w = 0.7 GeV and D = 1.4 
GeV^, for which the spectrum is plotted in Fig.[^ In all 
figures, our results are represented by blue circles, while 
experimental data are displayed as wide red boxes, where 
the height of the box indicates the size of the experimental 
error on the mass. Our error bars, where relevant, come 


from extrapolated results in situations where propagator 
singularities prohibit a direct calculation; details on the 
source of this problem and our extrapolation strategy can 
be found in the appendices of 

The next step is to identify a set of non-exotic splittings 
that provide similar results for our fitting-attempts, i. e., a 
reasonable correlation to the exotic case. In search of such 
a set we have found that axial-vector states seem to be 
good indicators. This is not entirely surprising; for exam¬ 
ple, as mentioned above, Ref. |45] focused on states from 
the J = 1 family and obtained a good overall description, 
including exotic states, using a separable kernel for the 
BSE, while both axialvector and exotic J =1 states had 
been difficult to describe otherwise. Another indication 
for the importance and connection to axialvector states 
is the hadronic decay of a tti into axialvector mesons and 
pions [21 1611162j , a feature which requires more work on 
our part and will be investigated in the future; still, it is 
important to note already here that a simultaneous good 
description of axialvector, exotic vector, and pion states 
is an excellent basis for a reliable subsequent calculation 
of their hadronic transition amplitude. 

With this in mind, we studied several sets of splittings 
and found that the combination of the ground-state pseu¬ 
doscalar to each of the axial-vector ground states together 
with the intra-axialvector splitting provide best guidance. 
The result is uj = 0.7 GeV and I? = 1.5 GeV^. As in¬ 
dicated by the proximity of the best-fit parameter sets, 
the results are indeed very close to each other and we 
therefore use this configuration to attempt predictions of 
1 ^ states for other quark masses. 
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FIG. 2. Charmonium spectra, fitted to axial-vector-related splittings. Calculated values (blue circles) are compared to 
experimental data (red boxes) from [541156tl^ . The horizontal dashed line marks the lowest open-flavor threshold. 


V. RESULTS AND DISCUSSION 

First, we proceed to charmonium, where the same fit¬ 
ting procedure yields uj = 0.6 GeV and D — 0.9 GeV^, 
for which the spectrum is plotted in Fig. In addition 
to our results (blue circles) and experimental data (red 
boxes), the lowest open-flavor threshold is plotted as a 
horizontal dashed line. With the pseudoscalar and ax¬ 
ial vector ground states well-met by the fit, we observe 
further adequate description of the scalar ground state. 
On the other hand, our calculation misses the vector and 
(pseudo-)tensor ground states, which is simply a conse¬ 
quence of the different anchoring and targets of the fit, 
as compared to our results in m- 

On a similar note, one observes that the level ordering 
of the ground states in the axialvector 1++ and 1'' , while 

correct for the light-quark case, is reversed in charmonium. 
This is a result of a rather straight-forward dependence 
of this particular splitting on the model parameter D: 
small values of D result in the correct level ordering, 
while increasing D at some point flips the level ordering 
and increases the splitting in the wrong direction. In 
our study, the requirement to reproduce the intra-axial- 
vector ground-state splitting is combined with a second 
requirement, which together yield a good description of 
the exotic light meson spectrum, namely to reproduce 
the pseudoscalar-axialvector splittings. The value of D 
found via the fit to this combination thus determines the 
axial-vector level ordering and splitting and can miss it 
by a small amount. Apparently, the charmonium result 
for D = 0.9 GeV^ is beyond the flipping value with the 
splitting still small, albeit reversed, which is acceptable 
for our purposes. 


The lowest charmonium state in the 1 channel has 
a mass of about 3.52 GeV, the next comes out at about 
3.88 GeV. 

The situation is similar in the bottomonium case, where 
our fit arrives at w = 0.7 GeV and U = 0.8 GeV^, for 
which the spectrum is plotted in Fig. While in this 
case our description of all measured ground states is 
rather accurate, our results for the excited states in many 
channels seem to be too low in general. Our resulting 
1 ^ ground state appears at 9.82 GeV and the next at 
10.02 GeV. 

While we set out from the isovector tti states, we now 
return to the light-quark domain for a brief comment on 
the isoscalar case. In our RL study, the BSE kernel does 
not account for flavor-mixing. While one can always mix 
the RL-BSE results afterwards EH, we refrain from doing 
so here, since there is no straight-forward and consistent 
way to implement such a mechanism for our present pur¬ 
poses. In a world with ideal SU (3) flavor mixing and 
assumed isospin-symmetry, an RL study provides the 
same results in the light-meson isovector and isoscalar 
channels. As a consequence, in such a world there should 
also be isoscalar versions of the tti states, based on our 
simple setup. In addition, at a slightly higher mass, there 
would be analogous states with ss quark content. 

As far as the interpretation of all of these new states is 
concerned, they may be hybrids or mix with other types 
of states, much like in other approaches to this feature 
of QCD. If one makes use of a hybrid interpretation 
of our results, the immediate question arises, how one 
would find the partners in possible corresponding hybrid 
supermultiplets. However, a study with the ability to 
identify hybrid and conventional states among our results 
is beyond the scope of our present investigation. On 
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FIG. 3. Bottomonium spectra, fitted to axial-vector-related splittings. Calculated values (blue circles) are compared to 
experimental data (red boxes) from [541156] . The lowest open-flavor threshold is above all masses plotted in this figure. 


another note, we would like to emphasize that our results 
for the 1 ^ channel, if compared to the results for hybrids 
presently available in the literature [5], are at the lower 
end of the apparent mass range: the predicted hybrid 1 ^ 
masses from lattice QCD, e. g., are substatially higher 
[T^ . whereas our results are in good agreement with, e. g., 
those from QCD sum rules [63) . 

While we have not aimed to describe states with the 
exotic quantum numbers = 0 , 0“' , or 2~^ , we 

do obtain results in these channels and we have included 
our results for these states in our figures for the sake 
of completeness. Considering the strategy in our setup 
specifically anchored to the J = 1 channels, we have 
no immediate reason to expect that our results present 
reasonable predictions. It is well conceivable that for the 
various exotic sets of quantum numbers there are different 
mechanisms at work and different parts of the interaction 
responsible for a good description. Naturally, one can 
turn to other approaches for guidance, which we may 
explore elsewhere. At present, our focus remains on a 
straight-forward study anchored to experimental data. 

A final word of discussion concerns the fact that the 
states in our setup are bound states and not resonances, 
since a decay mechanism is not incorporated in the BSE 
kernel in RL truncation. In particular above the rele¬ 
vant decay thresholds, this feature must be met with 
caution. In a first step, a semi-perturbative way to cal¬ 
culate hadronic decay widths is available and has been 
tested successfully in the past, see [64] and references 
therein. Since, in particular, the tti states are resonances, 
it remains to be seen whether such a treatment, using 
our current results, provides also hadronic widths of the 
correct size. Furthermore, considerable effects may also 


be expected from final-state interactions |65j . 

VI. CONCLUSIONS AND OUTLOOK 

We presented a covariant study of the masses of prospec¬ 
tive = 1 exotic mesons in both heavy and light 
quarkonia based on the DSBSE approach to QCD. Start¬ 
ing from a good description of the two lowest-lying tti 
states and identifying a good “marker”, namely the ax¬ 
ial vector channels, we transfer our assumptions to the 
other quarkonium systems and provide predictions for the 
masses of analogous states there. We see the lowest-lying 
1 state in both charmonium and bottomonium appear 
at about the same mass as the axial vector ground states. 
In summary, we have provided the hypothesis that, if our 
model assumptions are sufficient to capture the necessary 
physics and if the two lowest-lying tti states are well- 
represented by our setup, then one should expect to find 
analogous states also in the other quarkonium systems as 
described in detail above. 

To broaden our set of results we will investigate the 
decay properties of such states in the near future along 
the lines of previous studies of electromagnetic, leptonic, 
and hadronic transitions in our approach, which have 
been proven feasible and successful phenomenologically, 

e.g., in [saiMlISe]. 
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